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Field-Induced Gaps in the Frustrated Spin Ladder 
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We study the magnetization process of the 5 = 1/2 antiferromagnetic spin ladder in the 
presence of the second and the third-neighbor couplings which lead to frustration with the 
typical nearest-neighbor coupling. We use degenerate perturbation theory and level spectroscopy 
analysis of the numerical diagonalization data of the Hamiltonian for finite systems. We find two 
kinds of plateaux at half the saturation moment in the magnetization curve. One is mainly due 
to the second-neighbor couplings and the other to the third-neighbor couplings. The mechanisms 
of these two plateaux are quite different with each other. 
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The spin gap is a current topic of interest in strongly 
correlated electron systems because it is related to vari- 
ous interesting quantum phenomena such as high-T c su- 
perconductivity. Since the Lieb-Schultz-Mattis (LSM) 
theorem LP was recently generalized to the magnetization 
processEP , a field- induced spin gap has attracted a great 
deal of interest in the field of low-dimensional magnets. 
The extended LSM theorem predicts that a ID quan- 
tum spin system possibly has gaps which are observed 
as plateaux in the magnetization curve under the condi- 
tion of quantization of the magnetization, described as 



Q(S — to) = integer, 



(1) 



where Q is the spatial period of the ground state mea- 
sured by the unit cell. S and m are the total spin and 
the magnetization per unit cell, respectively. Applying 
this theorem to the spin ladder, only the well-known 
spin gapooB 1 is expected to appear at m = 0, as far as 
Q = 1. Several theoretical analyses, however, predicted 
that field-induced gaps would also appear at a finite mag- 
netization, with some modifications in the structure of 
the unit cell such as three-legtf and bond-alternating lad- 
ders.LP 

On the other hand, spontaneous breaking of the trans- 
lational symmetry (Q > 2) can also yield magnetization 
plateaux. The previous size scaling study© based on 
conformal field theory indicated the possibility of the 
plateau at m = 1/2 due to the two- fold degeneracy of 
the ground state (i.e., Q = 2) in the standard spin lad- 
der with the second-neighbor (2-N) interaction J2. In the 
plateau phase, the singlet pair state and the | ft) state 
of the rung are expected to locate alternately along the 
leg. A similar mechanism of the field-induced gap was 
predicted for the zigzag ladder equivalent -Up-i-hc bond- 
alternating chain with the 2-N interaction.li3'BE!P 

In this paper, we investigate another mechanism of the 
plateau formation due to the introduction of the third- 
neighbor (3-N) coupling J3. We also present a typical 



phase diagram of the J3-J2 plane, obtained by the level 
spectroscopy method analyzing the finite cluster diago- 
nalization data. 

We consider the 5=1/2 antiferromagnetic spin ladder 
with 2-N and 3-N exchange interactions in a magnetic 
field described by the Heisenberg Hamiltonian 



H = Hq + Hz 

L 

Hq = J\ 5^(S , l ) t ■ <S" + S2.1 ■ S24+1) 

i 

L 

+J± £ Sl.i ■ S2,i 

i 
L 



(2) 



Sl,i+2 + S2 



, 2,i+2) 



(3) 



(4) 



under the periodic boundary condition, where Ji, J±, 
J2 and J3 denote the coupling constants of the leg, rung 
and 2-N (diagonal) and 3-N exchange interactions, re- 
spectively (Fig. 1). Hereafter we put J±=l. Hz is the 
Zeeman term where H denotes the magnetic field along 
the z-axis and the eigenvalue M of the conserved quan- 
tity Si i + ^ii) i s a good quantum number. The 
macroscopic magnetization is represented by to = M/L. 

In order to consider the possibility and the mechanism 
of the magnetization plateau at m — 1/2 we use the 
degenerate perturbation theoxv around the strong rung 
coupling limit J\, J2, J3 <C IB We introduce a pseudo 
spin T for each rung pair and map the two original states 
(I U) - I lt))/\/2 and I H) of the S picture to the | 4) 
and I •ft") states of T, neglecting the other two states 
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Fig. 1. Spin ladder with 2-N and 3-N exchange interactions along 
the diagonals. 



(I T4> + I 4T))/V2 and | J.J.). After the mapping, we 
obtain the effective Hamiltonian 

L 



T L 



rpX 



T, 



i + 2) 



(5) 



This effective Hamiltonian describes the T = 1/2 XX Z 
chain with 2-N interactions, where Ji and J3 control the 
XX Z anisotropics and 2-N couplings, respectively, with 
fixed J\. We note that the XX Z anisotropy parame- 
ters of the NN and 2-N interactions are different from 
each other. Well-established woxks on this model have 
revealed the following properties x3 the system has three 
phases; the spin fluid (gapless), Neel (gapful) and dimer 
(gapful) phases. The J 2 = J3 — case is clearly in 
the gapless phase and sufficiently large J 2 {J3) yields 
the Neel (dimer) phase via, .the Berezinski-Kosterlitz- 
Thouless (BKT) transition.li-jlia) The boundary between 
the Neel and dimer phases is the Gaussian line. The 
above properties lead us to the conclusion that J2 and 
J3 give rise to the plateau at m — 1/2 in the original sys- 
tem, based on different mechanisms; the Neel state and 
the dimer state, in the language of pseudo spins. The two 
plateaux are hereafter called the plateau A (Neel) and 
plateau B (dimer). Clearly, both plateau phases should 
be accompanied by the two-fold degeneracy due to the 
spontaneous breaking of the translational symmetry. 

Next, we perform a numerical analysis for the original 
Hamiltonian (4) to more quantitatively confirm the real- 
ization of two plateaux at m — 1/2, predicted by the de- 
generate perturbation theory. The level spectroscopy^ 
is a powerful method to determine the BKT boundary, 
as well as most second-order transitions in ID quantum 
systems. In this method, the phase boundaries can be 
determined from the level crosses between low-lying ex- 
citations. This method is free from the most dominant 
logarithmic size corrections, which make it difficult to de- 
termine the BKT boundaries when conventional methods 
are applied. Hereafter, E(L,M,k) indicates the lowest 
eigenvalue of the Hamiltonian Hq in the subspace where 



the eigenvalue of J2i i^ti + is M and the momen- 
tum is k for the system size L. Using the Lanczos al- 
gorithm, E(L,M,k) is calculated for L = An (< 16), to 
avoid frustration among the 3-N exchange interactions 
under the periodic boundary condition. Before inves- 
tigating the full Hamiltonian, let us briefly demonstrate 
the powerfulness of level spectroscopy by drawing a phase 
diagram for the J3 = case, which has already been ob- 
tained by conventional methods.EP In level spectroscopy, 
we use the following two excitations given by 



Ai 



E L 



L 



1,7T 



E L 
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1,7T 



and 



(0) 



(7) 



When the parameters are swept, the state is gapless or 
gapful (i.e., plateau) according to whether_Ai < A or 
Ai > Ao, as explained by Okamoto et alS3 The phase 
diagram of the Ji-J\ plane, obtained by this procedure 
is shown in Fig. 2. The estimated error bars of the 
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Fig. 2. Phase diagram on the J2-J1 plane at m = 1/2. 



boundary points are much smaller than the size of the 
marks. In this phase diagram, in addition to the original 
spin fluid (gapless 1) phase, there appears another gap- 
less phase (gapless 2). Here we do not touch the gapless 
2 phase (equivalent to the S = 1 chain at m = 1/2; see 
refs. 8 and 17 for details). Using the effective Hamilto- 
nian eq. (5), we see that the slope of the boundary line 
between the plateau A and gapless 1 phases is 1 /3 in the 
limit of Ji — > (note thaL_f/eff of eq. (5) is exact in this 
limit), as noted by Mila.E3. In our numerical calcula- 
tion, we obtain J2/J1 — 0.3350 when J\ = 0.01, which 
agrees very well with the exact value 1/3. This shows 
the high reliability of our level spectroscopy method. 
Whoa. jConventional methods are applied to this prob- 
lem,E-3'B ) the slope value in the J\ — > limit is estimated 
as J2/J1 — 0.45, which is much larger than the exact 
value l/3—.The reason for this difference is discussed by 
OkamotoES* in detail and will be published elsewhere. 
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Here wc consider the full Hamiltonian problem J3 7^ 0. 
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Fig. 3. Excitations AqAi Aqb and Ai. 



Fig. 4. Phase diagram of the J2-J3 plane with fixed J\ (=0.4) at 
m = 1/2. 



For simplicity, we fix the 2-N interaction as J\ = 0.4. 
The gapless 2 phase does not appear in this case. Here, 
we consider two excitations having M = L/2 and k = n, 
because the change in the lowest energy level in this 
sector occurs when J3 is increased. Then, we define 
Ea(L, ■§, 7r) (Eb{L, -§, 7r)) as the lower energy in the re- 
gion of large J 2 (J3) and small J3 (J2), and also define 
A a (Aqb) in the same way as the eq. (7). These two 
excitations Ea(L, 4,7t) and Eb(L, 4,7r) correspond to 
the Neel and dimer excitations (see ref. 13) in the pic- 
ture of the pseudo spin T, and can be distinguished by 
the eigenvalues P of the space inversion operation of the 
rung number i — > L — i + 1. Namely, when L = An, 
the state for Ea(L,^,it) has P = — 1 and that for 
Eb(L, j,tt) has P = +1. In the spin fluid (gapless 1), 
plateau A and plateau B phases, the lowest excitations 
should be Ai, Aoa and Aob, respectively. Therefore, 
the phase boundaries can be determined frorn_the level 
crossing points among these three excitations.EiJ Figure 
3 shows these three excitations as functions of J3 when 
J\ = 0.4, J 2 = 0.14, L = 12. Thus, the state is gapless 
for J 3 < 0.040, plateau A for 0.040 < J 3 < 0.095, and 
plateau B for J 3 > 0.095. Repeating such procedures 
with sweeping the parameters, we obtain the phase dia- 
gram of the Ji — J3 plane as shown in Fig. 4. In order to 
investigate the universality of the boundary between the 
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— (L — > 00). (8) which strongly confirms the universality class of 



Fig. 5. Central charge c and critical exponent r\ z near the BKT 
boundary between the gapless and plateau A phases, as well as 
excitations Ai and Aqa- 



roles of Aoa and Aob are interchanged near the gapless 
1-plateau B boundary. Since the most dominant loa 
rithmic size corrections are canceled out in eq. ( |l0|),li 
we can obtain an accurate value of r\ z from eq. QlOj) . 
At the BKT transition point of the present type, the 
exponent rj z should be unity. Figure 5 shows the behav- 
iors of c and rj z near the gapless 1-plateau A boundary, 
as well as the level cross between Ai and Aoa, result- 
0.15. We can clearly see that r/ z ~ 1 at 
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where v s is the sound velocity, which can be calculated 

by 
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,(L^oo)(9) 



We also estimate the critical exponent rj z defined by 
(S$S?) - (S z ) 2 - (-l) r r-"* by use oM 



A 0B ) (10) 
near the gapless 1-plateau A boundary. We note that the 



^ = 2^(3A 0A 



the BKT transition. The behavior of c also suggests the 
BKT transition. A similar conclusion is also obtained 
for the gapless 1-plateau B transition. 

To clarify the properties of the boundary between 
plateaux A and B, we show the J3 dependence of the 
scaled gap 2LAi (2Ai is the length of the plateau of 
system size L) along the line J2 = 0.3 in Fig. 4. The size 
dependence of the scaled gap suggests that the plateau 
is opening in both phases and that the system is gapless 
only at the boundary J3 ~ 0.15. We also found c = 1 on 
the line labeled by open squares in Fig. 4. These results 
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In summary, the magnetization process of the frus- 
trated spin ladder was investigated with degenerate per- 
turbation theory and the level spectroscopy. The present 
analysis revealed the appearance of a novel magnetiza- 
tion plateau at m — 1/2 due to the 3-N interaction. The 
mechanism of this plateau is explained by the sponta- 
neous dimerization of the pseudo spin system. As far 
as we know, this is the first theoretical finding of the 
change in the plateau mechanisms, both of which have 
the spontaneous symmetry breaking, by sweeping phys- 
ically natural parameters. 
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Fig. 6. Behavior of the scaled gap 2LAi in the Ji=0.4 and J2 = 
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are consistent with the Gaussian fixed line predicted by 
degenerate perturbation theory. 

Several magnetization curves are also presented on line 
J 2 = 0.3; J 3 =0, 0.10, 0.15 and 0.20. They are obtained 
by the size scaling analysisEJ basied on conformal field 
theory and Shanks transformationO , using the numeri- 
cal result of E(L, M, k) up to L = 16. Only the lines of 
fitting polynomials are shown in Fig. 7. This suggests 




Fig. 7. Magnetization curves for Ji = 0.4 and J2=0.3, with var- 
ious J 3 (=0, 0.05, 0.1, 0.15 and 0.2). 
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that, with increasing J3, the plateau decreases until van- 
ishing at the critical point and then increases again. This 
behavior also explains why the mechanism of the gap for- 
mation due to J3 is different from the one due to J2. 

According to the present analysis, the gapless-plateau 
critical value of J3 for J2 — is smaller than that of J2 
for J3 = 0, irrespective of J\. In addition, the plateau 
B phase can appear for any ratio of J\/Jx> although the 
plateau A phase cannot appear for J±/J± > 1. Thus, the 
J3-induced plateau is more realistic than the pae due 
to J 2 . In fact, a typical spin ladder, SrC^OsEP, was 
reported to hold J\jJ\_ ~ 2. Thus, a plateau might 
be caused by the 3-N coupling of J3 in this or related 
materials. 



